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Scalar modifications of gravity have an impact on the growth of structure. Baryon and Cold 
Dark Matter (CDM) perturbations grow anomalously for scales within the Compton wavelength 
of the scalar field. In the late time Universe when reionisation occurs, the spectrum of the 21cm 
brightness temperature is thus affected. We study this effect for chameleon-f(R) models, dilatons 
and symmetrons. Although the f(R) models are more tightly constrained by solar system bounds, 
C^l " and effects on dilaton models are negligible, we find that symmetrons where the phase transition 

occurs before z* ~ 12 will be detectable for a scalar field range as low as 5 kpc. For all these models, 
the detection prospects of modified gravity effects are higher when considering modes parallel to the 
line of sight where very small scales can be probed. The study of the 21 cm spectrum thus offers 
a complementary approach to testing modified gravity with large scale structure surveys. Short 
scales, which would be highly non-linear in the very late time Universe when structure forms and 
where modified gravity effects are screened, appear in the linear spectrum of 21 cm physics, hence 
CN| 1 deviating from General Relativity in a maximal way. 
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■ I. INTRODUCTION 

A major challenge for theoretical cosmology is the explanation of the recent acceleration of the Universe's expan- 
+3 . sion [1]. In the standard A-CDM scenario, it is the consequence of the existence of a cosmological constant, although 
it could be also due to a dark energy fluid whose origin has yet to be determined [2]. Models of modified gravity [3] 
complement dark energy scenarios and provide an explanation of the absence of long range fifth force effects in the 
solar system and laboratory experiments. Indeed most involve at least one scalar field coupled to matter, and eventu- 
ally an environmental dependence leading to a screening mechanism of the scalar field in high density regions [4] . This 
mechanism is an essential ingredient for the models to pass the stringent constraints on the possible modifications of 
gravity in the laboratory [5], the solar system [6], and the galactic environments [7]. Moreover, the scalar fields are 
£s1 ■ required to sit at the minimum of the density dependent effective potential prior to Big Bang Nucleosynthesis (BBN) , 
so that catastrophic modifications in the formation of light elements are avoided. 

Numerous models of this type have been proposed. Let us mention the chameleons [8-15], involving a thin shell 
shielding the scalar field in dense bodies, the symmetrons [16-22], involving a symmetry breaking potential so that 
I the scalar field is decoupled from matter at high densities, the dilatons [23, 24], where the coupling to gravity turns 
off in dense environments, and the f(R) models [24-34], which are a sub class of chameleon models [24]. 

At the homogeneous level, all these scenarios coincide with the A-CDM model. However, the evolution of linear 
perturbations differ. As a consequence, models of modified gravity can induce observable signatures in the matter 
1 power spectrum at redshifts z < 2, and less importantly in the cosmic microwave background [24, 35] at z ~ 1100. 
Modified gravity models can also be probed with weak lensing (see e.g. [36]). During the dark ages and the reionisation 
. . . ■ period, i.e. in the range 1100 > z > 6, no cosmological signal have yet been observed. Such observations would be 
however of great interest for cosmology, especially for the study of modified gravity through the time evolution of the 
matter perturbations. 

In the near future, this gap is expected to be partially filled with the observation of the 21cm signal from reioni- 
sation [37-41], and maybe in the more distant future, from the dark ages. During reionisation, transitions between 
the fundamental hyperfine levels of neutral hydrogen atoms are possible, via the Wouthuysen-Field effect involving 
the absorption and re-emission of Lyman-a photons from the first stars (for a review, see Ref. [40]). These induce 
the 21cm signal corresponding to a stimulated emission of 21cm photons against the Cosmic Microwave Background 



> 
CO 



* Email address: 
t Email address: 
t Email address: 



philippe.brax@cea.fr 

s . clesse@damt p . cam .ac.uk 

a.c.davis@damtp.cam.ac.uk 



2 



(CMB) radiation. The 3D power spectrum of the 21cm radiation maps the baryon distribution and thus is sensitive to 
modifications of gravity. In this paper, we explore for the first time the effects of modified gravity on the 21cm power 
spectrum at reionisation, and discuss the detectability of such effects with instruments of future generation giant Fast 
Fourier Transform radio-telescopes [42, 43]. 

Several parameterisations of modified gravity have been proposed [44-56] for the evolution of linear perturbations 
such as, for instance, through the Poisson equation, — fc 2 $ = 411(1 + v)Ga 2 5p m and ^> = (1 + 7)$, where Sp m is the 
matter density perturbation and where $ and "J are the potentials in the Newtonian gauge. This parametrization 
involves two functions v(k, a) and j(k, a) that depend both on time and on the perturbation wavenumber k. For f(R) 
models, those functions depend on B = (/rb/ fn)HdR/dH together with fn today [32]. In this paper, we adopt the 
parameterisation proposed in Refs. [57, 58], for which the action and dynamics can be fully and uniquely reconstructed 
from two time-dependent functions: the coupling to matter /3(a) and the scalar field mass m(a). Of course, in the 
f(R) case, it coincides with the usual approach. Moreover, for all these models, the v and 7 functions can be explicitly 
obtained as a function of m(a) and /3(a). 

Using the (m(a),/3(a)) parameterisation, we calculate the signatures of f(R), dilaton, chameleon and symmetron 
models on the 21cm power spectrum at reionisation, as well as on the present matter power spectrum in the linear 
approximation. For each model, we discuss the range of parameters that could be probed by 21cm experiments and 
compare it to the constraints from local experiments. Finally, we evaluate how the coupling to photons could be 
bounded via 21-cm constraints on the variation of the fine structure constant a. In all cases we find that the 21 cm 
signal obtained by varying the modes parallel to the line of sight is the most relevant. In the f(R) case, we expect 
the constraints to be less stringent than the ones from the gravity tests in the solar system. For dilatons, the signal is 
found to be negligible while symmetron models with a transition at a rcdshift larger than ~ 12 could be detected 
even when the range of the symmetron interaction now is as low as 5 kpc. 

This paper is organised as follows: In Sec. II, we briefly summarise the physics of the 21cm signal at the period of 
reionisation and derive its 3D power spectrum. We also give the evolution of the baryon and dark matter perturbations 
after the time of last scattering for the A-CDM model. In Sec. Ill, we show how these equations are modified in the 
context of scalar models and describe the dynamics of f(R), symmetron, chameleon and dilaton models using the 
reconstruction of the scalar field dynamics from the parametrization of Ref [57, 58]. In Sec. IV we give the specifications 
of the considered FFTT radio-telescope as well as the forecast errors on the 21cm power spectrum for single redshift 
measurements. We then evaluate the range of parameter values leading to observable effects and compare to the 
constraints from local tests of gravity. Our results are summarised in the conclusion. 



II. 21CM SIGNAL FROM REIONISATION 



In this section, we briefly review the basics of the 21cm cosmic background from the period of reionisation and 
refer to Refs. [40, 41, 59, 60] for a more exhaustive description of the signal and its relevance to cosmology. The 21cm 
signal corresponds to an absorption (or a stimulated emission) of 21cm photons from (against) the CMB, induced 
by the transitions between the hyperfinc ground state of neutral hydrogen (HI) atoms. During the dark ages, after 
the thermal decoupling of the baryon gas at z ~ 200, the hyperfine level population is shifted away from thermal 
equilibrium with photons because of the spin changing collisions between HI atoms, and between HI and free electrons. 
In this paper, we focus on the 21cm signal from the period of reionisation [z ~ 10). During this period, collisions 
are rarefied and hyperfine transitions are driven by the so-called Wouthuyscn-Ficld effect, via the absorption and the 
re-emission of Lyman-a photons coming from the first luminous objects that reionise the Universe. 



A. 21cm Power Spectrum 



At the homogeneous level, the difference between the observed brightness temperature and the CMB temperature 
at a given observed energy E is given by [59] 

T B (E) = (l-e^')(^^ , (1) 

where 

3c 3 n m (a E )Ai h p a E 
TE ~ 327rr s (a B )fcB^ 2 1 (d^||/dr) ' [ ' 

is the 21cm optical depth. i/ 2 i is the frequency of the 21 cm line in the rest frame, = 1/(1 + Ze) is the scale factor 
for which E = a^E^i with E21 = h p U2i, the energy of the 21 cm spin flip transition. Aiq = 27r<w| 1 /ip/(3c 4 m 2 ) ~ 
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2.869 x 10~ 15 s" 1 is the Einstein coefficient of spontaneous emission, dv«/dr is the gradient of the physical velocity 
along the line of sight and r is the comoving distance. At the homogeneous level, there is no peculiar velocity and 
thus dvn/dr — cleH^e)- «hi is the number density of neutral hydrogen atoms. It is related to the baryon number 
density rib and the neutral fraction ih via nn/ = £H"-t>(l — /hc), where /n e is the Helium fraction. 

During the reionisation process, te *C 1, so that the 21 cm brightness temperature is well approximated by 



. , 3c 3 XH(oB)n b (a£;)(l - /hc)-4io^ p T s - T 7 
1b(&) = 



mvkBV^H \a E ) T s (l + z) 



(3) 



After perturbing at the linear level the baryon density, the ionized fraction and the gradient of the radial velocity, 
the brightness temperature in a given direction e is given by 

T B (E, e ) = ^ [1 - Xi(l + Si)} (1 + S h ) (1 - 5y) , (4) 

where a bar denotes homogeneous quantities. <5b and Si are respectively the relative baryon and ionized fraction 
perturbations, and 

S v = —-l 5 
all or 

is the radial gradient of baryon peculiar velocities. Let us note that in order to obtain Eq. 4, we have followed Rcf. [60] 
and assumed that there exists a redshift range during the reionisation with T s T 1 so that (T s — T 7 )/T s ~ 1. With this 
assumption, the fluctuations of the spin and photon temperatures are typically second order and can be conveniently 
neglected. 

Defining by 9 = diV % the divergence of the velocity field we have 

dv r ,q\ 
dr k ' 

where n is the line of sight vector and k the wave vector of a linear perturbation. The perturbed conservation of 
matter gives 



V, 



= -K (7) 



in conformal time. In the long time regime, the baryon contrast grows with a growing mode D+ which is not the scale 
factor a anymore in modified gravity, hence 

S b = D+S b0 , (8) 

from which we get 

drj 



dlnD, 



= S b , (9) 



and therefore 



1 /dlnD+\ 

^=-^(^rJ^' (10) 

where jl = k • n/fc is the cosine of the angle between the line of sight and the wave vector k. The first factor 
(din D '_)_) / ' {aH drf) is equal to one in the absence of modified gravity. In the following we will use 

H^r)(") 3 - 

From Eq. (4), the brightness temperature perturbation in Fourier space is given by 

AT B (k) = ^ {(1 + M 2 )«5 b - Si ^ - (1 + M 2 )4] } • (12) 

XH 
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Therefore the 21cm 3D-powcr spectrum reads [60] 

PatM = (J^J { [x^P hh {k) - 2x H x i P ih (k) + zffli(fc)] + 2/i 2 [x\P hb {k) - 5 H 5iPib(fc)] + y, A xlP hh {k)} , (13) 

where the power spectra P a p{k) with a,/3 = i,b arc defined from (<5 Q (k), 8p(k')) = (27r) 3 i5 3 (k — k')P a p(k). 

From Eq. (13), one can notice that the fj, 4 component does not depend on P^ and Pa. By measuring the 3D power 
spectrum, it is therefore in principle possible to separate cosmology from the astrophysical contaminants related to 
the reionisation process. In this paper, we will assume for simplicity an optimistic reionisation model [60] for which 
Pib = Pi = at the redshift of interest. 

A dedicated 21cm experiment does not directly measure the comoving modes k but the angular positions in the sky 
and the signal frequency. It is therefore convenient to determine the observable 21cm power spectrum in the Fourier 
dual of the space of angular positions and frequencies, the so-called u-space. It is related to the k-space through the 
relations 

u ± = D A k ± , (14) 

«|| = V k \\ , (is) 
where £>a is the comoving angular distance, given in a flat Universe by 

D ^ = c I ^ dz ' ^ 

and where y{z) is the conversion factor between comoving distances and frequency intervals, 
The 3D power spectrum of the 21cm brightness temperature in u-space reads 

< 18 > 

B. Evolution of Baryon Perturbations 

When baryon and dark matter perturbations re-enter inside the horizon, they feel their combined gravity. Baryon 
perturbations feel also the pressure. Assuming a A-CDM model at the background level, cold dark matter and baryon 
perturbations evolve according to [61] in the absence of modified gravity 

s c + 2hs c = hi 2 (n h s h + n c 6 c ) , (19) 

5 h + 2H5 h = ' A -H 2 (n b 6 b + n c S c ) - ^l^Is ^ + ^ ^ (20) 
2 cr fx m 

where f2b is the baryonic fraction which is a constant / = Pb/(Pc + Pb) in the matter era, f2 c the CDM fraction, /x m is 
the mean molecular weight and St is the gas temperature perturbation. The homogeneous gas temperature T g evolves 
during the dark ages according to [61] 



where 

3m 



t- 1 = — = 8.55 x 10~ 13 yr ■ (22) 
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The last term on the r.h.s. of Eq. (21) accounts for the energy injection due to the Compton scattering between CMB 
photons and the residual free electrons. During the period 1100 > z > 200, the Compton heating drives T g — » T 7 . 
After z ~ 200, due to the expansion, the gas temperature decouples from the radiation and evolves as T g oc 1/a 2 , as 
expected for an adiabatic non relativistic gas in expansion. We have assumed that the non-trivial evolution of the gas 
temperature during the reionisation has only a negligible effect on the baryon perturbations. 

Eq. (21) can be perturbed at the linear level. On small scales and assuming that photon density and photon 
temperature perturbations can be neglected, the gas temperature perturbations evolve according to [61] 



3 try d Tcy 



5t 



(23) 



At high redshifts, the interaction p + e <-> H + 7 maintains the species in equilibrium and the free electron fraction 
is given by the Saha equation 



1 



2tt 



3/2 



-eo/T 



(24) 



where eg = m c 



mjj = 13.6 eV. Near the redshift of the last scattering, the equilibrium is not maintained 



anymore and the Saha equation becomes inaccurate. One therefore needs to solve the Boltzmann equation for x e . A 
good approximation 1 is given in Refs. [63, 64], 



dx c 



(1 - x e )j8i - x\n\>a 



(2) 



(25) 



where 



2tt 



3/2 



-eo/T 



(26) 



is the ionization rate, and where 



v(2) 



9.78 



a- 



(?)" 2 '» 



(?) 



(27) 



is the recombination rate. The superscript ^ indicates that recombination in the ground state is not relevant. Indeed, 
this process leads to the production of a photon that ionizes immediately another neutral atom and thus there is no 
net effect. For the purpose of this work, we have calculated the free electron fraction with Eq. (24) for x c > 0.99, and 
with Eq. (25) for x c < 0.99. 

Equations (19), (20), (21), (25) and (23) form a closed set of equations. In order to calculate the 21cm power spec- 
trum, these have been integrated numerically. Initial conditions for the baryon, cold dark matter and gas temperature 
perturbations are provided at a redshift z = 950 by the CAMB code [65]. In the next section, the equivalent equations 
to Eqs. (19) and (20) for modified gravity models will be described. We have only considered scenarios for which 
modified gravity effects at z > 950 are negligible. 



III. MODIFIED GRAVITY MODELS 



A. Scalar Field and Modified Gravity 



Scalar-tensor theories are characterised by their coupling to matter and their interaction potential, 



S = J y^d 4 



(28) 



1 For more accurate results, the recombination can be calculated numerically by using the RECFAST code [62] . taking into account additional 
effects like Helium recombination and a 3-level atom. 
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where <p is the scalar field and V((f>) its potential, S m = S m g^u\ the matter action with ty 1 the matter fields which 
are minimally coupled to the Jordan frame metric g M „ = A 2 ((f))g fliy ; g^ v is the Einstein frame metric, which is used to 
compute the Ricci scalar R; n\ = k = 8nG = m~y where G is Newton's constant and m p i the reduced Planck mass. 
The field equations are obtained by varying the action S with respect to the field </>, and we have 

f = g^f^, (29) 

= jL* Sm , (30) 
v-9 og^v 

where we have defined the Jordan frame energy momentum tensor that is related to the Einstein frame one by 
T M „ = A 3 ((j))T fip gp I y. The Klein-Gordon equation for </> and the Einstein equations become 



n<f> = v g,4<i>,T) , (3i) 

(32) 
(33) 



nr _ \ Rg ^ = kT z , (32) 



in which the scalar field is governed by an effective potential 

V cS (^T) = V{4>)-A{4>)T , (34) 
and the total energy momentum is given by 

TZ = A(cf>)T^ - g^V(4>) + - (V0) 2 , (35) 

which satisfies the following conservation equation 

V^ v = ^ (Tg^ - T^) . (36) 

Note that this implies that for pressureless matter with T^ v = p m u^u v and u M the 4-velocity (u^u u = —1), we have 
— T = p m and p m is conserved independently of <p: 

V^u") = , (37) 

or equivalently 

p m + 3hp m = , (38) 
where p = u^V M/ o and 3 ft. = V^u^ 1 . The continuity equation is 

v? + Ki P(j)u^ = -niPd^cj) , (39) 
where the coupling to matter is defined to be: 

13 = Tripsin A . (40) 

We are interested in models of modified gravity with screening properties in dense environments due to the non- 
linearity of the interaction potential V and the coupling function A. Examples of such models are chameleons, dilatons 
and symmetrons. f(R) models are chameleons in disguise written in the Jordan frame. 



B. Perturbations 



We will focus on models where gravity is only modified at late times well after last scattering. The perturbation 
equations follow from the conservation of matter and continuity equations. After last scattering, baryons and photons 
decouple while the baryons and CDM start evolving in a coupled way under the influence of gravity. The conservation 
equations become 



(41) 
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and 

while the continuity equations are 

<9 b + W h = fc 2 $ + f3 h k 2 K 4 5(/) + k 2 c 2 S h , 
where we are taking into account the baryonic speed of sound Cb, and for CDM 

e' + ne c = k 2 ® + f3 c k 2 n±8<j) . 



(42) 
(43) 
(44) 



We have defined the divergence of the velocity fluid for both fluids 6 = d l Vi and generalised the models by considering 
that baryons and CDM couple differently to the scalar field. The Poisson equation is now 



~n 2 (n h 5 h + n c 6 c ) 



Similarly the scalar field satisfies the Klein-Gordon equation and in the sub-horizon limit becomes 



Ki8(j> = —iT-L z 



k 2 + m 2 a 2 



The baryon contrast satisfies 



elk 2 



\u 2 n h li 



5 h - ^n A n c i 



2/3 b /3c 



(45) 



(46) 



(47) 



and CDM 



H5' - -n 2 n h 1 



2/3 b /?c 



f 



s h - -n 2 n c i 



1 



4 = 



(48) 



These modified equations allow one to study the influence of modified gravity on 21 cm physics. The relative differences 
between the baryon perturbation evolution for the A-CDM model and for some models of modified gravity have been 
plotted in Fig. 1 for the wavelength mode k = O.lMpc" 1 in the observable range of 21cm experiments. 



C. Some Models 



1. Chameleons 



We will consider chameleon models where the coupling to matter f3 is constant and the potential V{<j>) is a decreasing 
function of <j). As a result the effective potential 

V en {4>) = V((f>) + e M ' m ^ Ph + e^' m ^ Pc , (49) 

where the couplings to baryons and CDM are taken to differ, has a minimum <j) m m- When the mass of the scalar 
field at this minimum m is large enough and greater than the Hubble rate H, the minimum is stable since before 
Big Bang Nucleosynthesis (BBN), guaranteeing a small variation of the masses of particles in the Einstein frame, and 
therefore a negligible modification of the formation of matter during BBN. This condition must be satisfied also to 
avoid important deviations from A-CDM at the homogeneous level. 

Typical examples of potentials correspond to inverse power laws of the Ratra-Peebles types. Local experiments such 
as cavity tests of gravity or the Lunar Ranging test of the equivalence principle impose stringent restrictions on the 
mass of the scalar field now. We will come back to these restrictions later. 
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2. f(R) models 

Viable f(R) models are nothing but chameleon models with a definite value of the coupling (3 = 1/V6, and with 
the potential given by 

V{4>) = ™pi , (50) 

where fn = df /dR. The mapping between R and <f> is given by 

f R = e - 2/3 ^ mpl . (51) 
Typically we will be interested in large curvature models R> R*, where 

I R \™ 

f(R) = R + R a + R 1 (^-J (52) 

and Rq plays the role of a cosmological constant and n > defines the asymptotic expansion of the f(R) function. 
For these models, the range of the scalar force is constrained by the requirement itiq/Hq > 10 3 which springs from 
the loose constraint that galaxies such as the Milky Way should have a thin shell. 



3. Symmetrons 

For the symmetron model, the interaction potential and the coupling function are simply chosen such that 

V{4>) = V - \/ sym 4> 2 + \\ S y^ , 

= ! + (53) 
in which Vb is a cosmological constant. The effective potential can be rewritten as 

V cS {4>) =\(j%2- /4m) ^ + \>**<P ■ (54) 
Hence, in the symmetron model, as long as p m is high enough, namely p m > p+ where 

P* = f4ym M2 , (55) 

the minimum of the effective potential is at the origin (</> = 0). In contrast, in vacuum, the symmetry is broken and 
the potential has two nonzero minima: 

& = . (56) 

When the matter density docs not vanish, the minimum of the effective potential depends on p m : 



0min(/Om) = <t>+\ 1 ~ — 9 (p* - p m ) , (57) 

V p* 

where 9(x) is the Heaviside function. The minimum is not zero only for p m < p* and converges to 0* when p m vanishes. 
As long as </> <C M the effective coupling to matter reads 

M « ^ , (58) 

leading to the absence of modification of gravity in dense environments where the field vanishes. Indeed, the variation 
of the effective coupling with the matter density is: 

Anm(Pm) =P*(l-^\6 (p* - Pm ) . (59) 
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The symmctron model is designed to induce modifications of gravity which could be tested in the near future, 
both gravitationally and cosmologically. Requiring that the energy density at which the curvature at the origin of the 
potential changes sign is close to the current critical energy density , we have the estimate 

AfV'ym ~ #0»4l • (60) 

Moreover, the modification to gravity is detectable only if it is comparable to (or bigger than) standard gravity, or 
equivalent ly the effective coupling /? ~ 0(1), which implies 

M ~ ' 

where we have defined Ai = m P1 /M 2 . These determine the vacuum mass 

m 2 (0 o ) = 2 Ms 2 ym ~ O Hi = 0{A 2 )Hl (61) 

and correspondingly set the interaction range of the symmctron to be ~ O [m _1 (<^>o)] • From the study of solar system 
tests, it can been found that M < 10 _3 mpi or equivalently 

A 2 > 10 6 . (62) 

This implies that 0*/M < 10 -3 . The range of the symmetron in the cosmological vacuum is given by 

fi£m S iO^o" 1 ~ 10 Mpc , (63) 
which corresponds to relevant scales for astrophysics. 



4- Dilatons 



Dilaton models arc akin to symmetrons in as much as the effective potential has a minimum whose origin, here, 
springs from the coupling to matter 

' - <M 2 , (64) 

4>*) ■ (65) 

On the contrary, the bare potential V(<j>) is assumed to be a smooth and non- vanishing function whose order of 
magnitude V(4>+) corresponds to the vacuum energy now. The dynamics can be analysed close to 4> = 0* where the 
field is attracted in high density environments. Local constraints on modifications of gravity impose that 

M > 10 6 , (66) 

in a similar fashion to the symmetron case, which implies here that ttiq/Hq > 10 3 again. 



which leads to a universal coupling 



D. Reconstructing the Dynamics 



All the models which have been presented in the previous section can be parameterised solely using the time evolution 
of the mass and coupling functions at the minimum of the effective potential. This is a particularly convenient way of 
defining models when one is interested in cosmological perturbations. The dynamics of the models of modified gravity 
of the chameleon, dilaton and symmetron types are completely determined by the minimum equation of the effective 
potential 



dV 



= -13- 



mpi 



(67) 
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FIG. 1: Relative differences between the linear evolution of baryon perturbations with k = 0.1 Mpc -1 for the A-CDM model and 
for f(R) and symmetron models. Plain curves are for the symmetron model. From top to bottom lines, the model parameters 
are respectively mo = 0.01 Mpc _1 ,/3o = 1,2* = 20 (red curve), mo = 0.1 Mpc _1 ,/?o = 1,2* = 20 (yellow curve), mo = 
0.1 Mpc -1 , Pa = 1, 2* = 14 (green cuvre) and mo = 0.1 Mpc -1 , Po = 0.5, z* = 20 (blue curve). Dashed curves are for the f(R) 
model, with mo = 10 -4 Mpc -1 (top red dashed curve) and mo = 10 -3 Mpc -1 (bottom blue dashed curve). 



where we have introduced p = 3H m^, 



p = p c A c n c + p h A h n h 



(68) 



and Ab, c ~ 1 due to BBN constraints. In fact, the knowledge of the time evolution of the mass to and the coupling /3 
is enough to determine the time evolution of the field and the potential completely. Using the minimum equation we 
deduce that the field evolves according to 



at 



m mpi 



(69) 



This is the time evolution of the scalar field at the background level since the instant when the field starts being at 
the minimum of the effective potential. This leads to the solution 



0(a) 



3 
mpi 



/3(«) 
am 2 (a) 



p(a)da + cj) c , 



(70) 



where <p c is the initial value of the scalar field. Similarly the minimum equation implies that the potential can be 
reconstructed MS M function of time 



V = V - 3 



P(a? P 2 
am 2 {a) m 2 ^ 



da, 



(71) 



where Vq is the initial value of the potential at a = a^. This defines the bare scalar field potential V((j>) parametrically 
when P(a) and m(a) are given. Let us come back to the chameleon, f(R), dilaton and symmetron models for which 
one can explicitly verify that this way of defining models can be used. 

1. Chameleon and f(R) models 

For these models and in the matter era, the coupling to matter j3 is constant. We will be interested in the models 
where 



m = toq a 



(72) 
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When r > 3 and /3 = 1/Vo, they correspond to the large curvature f(R) models with r = 3(n + 2)/2 in the matter 
era. When 3/2 < r < 3, the models are of the chameleon type with an inverse power law potential 



v{4>) 



A"+ 4 



(73) 



and n={2r- 6)/(2r - 3). 



2. Symmetron 



The symmetron models can be reconstructed using 

P(a)=(3^l-(^) 3 (74) 
for z < z+ and /3 = 0, z > 2*. Similarly we have 

m(a) = m+^J 1 — ^ — . 
The parameters (/3*,m*,2*) determine the model completely, 



(75) 



,2 



2/?*p* . Msym 

m* = V2/i sym , A sym = — — . (76) 



Finally we have 



m = ^ ■ (77) 



3. Dilatons 



For dilatons, the behaviour of the mass and coupling functions in the matter era close to 4>* is universal and defined 

by 

m = 77io a_2 ! (78) 

corresponding to a mass which is proportional to the Hubble rate. The coupling to matter is time dependent and 
follows the inverse matter density 

13 = Poa 3 , (79) 

with a coupling which increases as matter becomes sparser. This model can be seen as a generalisation of chameleon 
models where one can choose 

m = moa~ r , ft = /3oa~ s , (80) 
where r = 2 and s = —3 here. We will consider these (r, s) models in the following. 



IV. 21CM POWER SPECTRUM OF MODIFIED GRAVITY 



A. 21cm FFTT-type experiment 



The present and next generations of large radio-telescopes, like LOFAR [66, 67], MWA [68] and SKA [69, 70] are 
designed for the detection of the 21cm signal from the period of reionisation. However, their sensitivity is not expected 
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to be sufficient to constrain cosmology 2 through the observation of the 21cm power spectrum [60]. Nevertheless, 
M. Tcgmark and M. Zaldarriaga have recently proposed the concept of Fast Fourier Transform radio- Telescopes 
(FFTT) [42] whose potential ability to measure the 21cm power spectrum and to put strong constraints on the 
cosmological parameters have been demonstrated in Rcf . [60] . The FFTT is an all digital radio-telescope composed of 
a square grid of dipole antennas. The multifrequency images of half the sky are reconstructed from the data measured 
by each antenna after several fast Fourier transforms. The key advantage of the FFTT compared to traditional 
interferometric radio-telescopes is its cost, scaling as Na log 2 Na instead of N^, with Na the number of dipole 
antenna. 

In this paper, we adopt the FFTT as a template to study modified gravity effects on the 21cm power spectrum from 
reionisation. We refer to Refs.[42, 60] for the main specifications of the experiment that are reported in Tab. I. The 
aim of this section is to evaluate the errors on the 21cm power spectrum for the considered FFTT experiment. For 
simplicity, we assume ideal foreground removals and refer the interested reader in foregrounds and removal techniques 
to Refs. [71-73]. 

The noise spectrum P n for the FFTT is given by [42] 



P n(u) = f2 Z o , ^K)^(^ll) , (81) 
/covcr^max 1 'fovEo 



4^/ sky A 2 T 2 
2 D 2 Slf 

"over max " 

where T sys is the system temperature, D max is the length of a side of the FFTT, f2f ov = 2ir is the total field of view, 
/sky = f2fov/47r and t Q is the total observation time. / COVC r = Na(D / D mayi ) 2 is the fraction of the total area covered by 

the antenna. D is the minimum baseline between two antennas. P^ and Pn are the fouricr transform in the u-space 
of the gaussian angular and frequency window functions. These are introduced to take into account respectively the 
angular resolution and the frequency resolution of the telescope. They read 

PnM = ^ W± > (82) 
Pj(«ll) = e"H , (83) 

with 

a ± = , an = Au , (84) 

A nax V8hi2 11 

and where Av is the frequency bandwidth. Moreover, to avoid non-linear effects, we follow Rcf. [60] and assume a 
sharp cut off at fc max = 2Mpc _1 , so that for our specifications of the FFTT, we can safely consider pj = 1. 
It turns out that the error on the power spectrum for a given mode u reads 

«Wu) = - J= [PAT b (u) + P n (u x )] , (85) 



where 

iVc = ^|^ (86) 

"patch AU_l 

is the number of cells in u-space probed by the experiment. For the fiat sky approximation to be valid, it is required 
to divide the sky in small patches whose size f2 pa tch is typically smaller than lsr, so that Auj_ = 2n/ y/£l pa ,tch- 

Because the FFTT is a signal-dominated experiment in the range of wavelengths of interest for cosmology, the error 
on the power spectrum is mainly due to the cosmic variance. It is noticeable that compared to CMB experiments 
which only probe the last scattering surface, the induced errors on the cosmological model parameters can in principle 
be reduced by the observations of the 21cm 3D power spectrum at several redshifts. 



B. 21cm Spectra and Explicit Models 

In this section, we calculate the effects of modified gravity on the 21cm power spectrum at reionisation and on the 
matter power spectrum today. For each model, we give rough bounds on the model parameters to have observable 



2 Except for the SKA and for some optimistic scenarios [60] 
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Total size 


fmax = 1 km 


Min. baseline 


D= 1 m 


Number of antennas 


iV A = io 6 


Bandwidth 


Az; = 0.05 Mhz 


System temperature 


T sys = 400 K 


Observation time 


t = 1 year 


Angular resolution 


reB = \/D (Beam FWHM) 


Field of view 


fifov = 27T 



TABLE I: Specifications of the considered FFTT experiment. A is the redshifted wavelength of the 21cm signal. We assume a 
gaussian beam and use the FWHM convention [42] for the beam width, as well as ideal foreground removal. The FFTT covers 
half of the sky sphere. 



signatures with the FFTT. Then these bounds are compared to those obtained with the matter power spectrum and 
with local tests of gravity. At the end of the section, we show that in principle the coupling of the scalar field to 
photons can also be constrained with 21cm observations through the variation of the fine structure constant a. 



1. f(R) models 

The relative differences between the 21cm power spectra at z = 11 for the f(R) model and for the A-CDM model 
are plotted in Fig. 2 for wavelength modes orthogonal to the line of sight (/x = 0), for r = 3, /3 = l/y/6, and various 
values of the parameter to . The case r=3 corresponds to the limit case n — > of large curvature f(R) models. The 
expected relative errors on the power spectrum measurements for the FFTT experiment arc also added. 

On large scales (k < O.OlMpc^ 1 ), the terms in /3 2 /(l + m 2 a 2 /k 2 ) in Eqs. (47) and (48) tend to zero and the f(R) 
model behaves like the A-CDM model. The relative difference with the A-CDM model grows at small scales, where 
these terms become important. For sufficiently small values of too (< 10~ 4 Mpc _1 ), the 21cm power spectrum can 
be strongly affected by modified gravity effects. For these very low masses, the mass at z = 11 is still larger than 
the Hubble rate. This is not the case now implying that the scalar field does not follow the minimum of the effective 
potential in the recent past of the Universe, leading to potential deviations from A-CDM at the background level at 
small rcdshift. These cases are already excluded by local constraints and are only presented for illustration. A difference 
with the A-CDM could be detected by a FFTT experiment up to too ~2x 10~ 3 Mpc -1 for modes orthogonal to the 
line of sight. 

The situation is improved if we consider k± w O.lMpc -1 (in order to have a good compromise between a large value 
of N c and ~ 1), and k\\ up to its maximal value fc max = 2Mpc _1 (remember that the noise is independent of fcy). 
The last bound can be enhanced by one order of magnitude in this case, up to too ~ 2 x 10 _2 Mpc _ , as shown in 
Fig. 3. Hence, by using the frequency component of the 3D-21cm power spectrum, one probes smaller perturbation 
scales for which the effects of modified gravity arc potentially more important. Such modes cannot be probed directly 
with the matter power spectrum today since they correspond to non-linear scales for which the screening mechanism 
suppresses the deviations from the A-CDM behaviour. 

This limit on the scalar field mass corresponds to a ratio mo/Ho ~ 90, i.e. lower than the constraint obtained with 
local tests [57], mo/Ho > 10 3 (imposed to guarantee that galaxies such as the Milky Way have a thin shell). However, 
it is comparable to the solar system constraint [57], 

to| /3 2 a n olO-^+ 12 

(2r-s-3)$ ' {S<) 

where <f> ~ 10~ 6 is the solar Newtonian potential, giving rise for s = 0, r = 3, j3 = l/y/6 to the bound m /H > 100. 

We have also calculated the matter power spectrum today for too — 2 x 10 _3 Mpc~ 1 in the linear approximation. It 
is plotted in Fig. 8. The matter power spectrum is observed to be slightly outside the 68% C.L. error bars of the SDSS 
data. For increasing values of the scalar field mass, the difference with the A-CDM matter power spectrum becomes 
undetectable. However, it must be noticed that we did not take into account potentially important non-linear effects. 

To summarise, we have found that future observations of the 21cm signal at reionisation with a FFTT-like radiote- 
lescopc should constrain more efficiently the f(R) model with r — 3 than the present matter power spectrum, mainly 
because the signal can probe smaller scales in the linear regime. Modified gravity effects could be detectable up to a 
scalar field mass value mo ~ 2 x 10 -2 Mpc -1 , corresponding to ttlq/Hq ~ 90. Although this value is lower than the 
galaxy constraint from [57], it is competitive with the tests of gravity in the solar system and better than the bounds 
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from the CMB and other local tests of gravity. It will therefore be interesting to study whether this situation could be 
improved by using multi-redshifts measurements and by combining data from the whole 3D u-space, or made worse 
due to possible degeneracies with cosmological and nuisance reionisation parameters. Such a study will require ideally 
the use of more complex Fisher matrix or Monte-Carlo methods and is left for future work. This last remark is also 
valid for the other models we consider in this paper. 




-0.2 ■ ■ ■ — . . ■ . i ■ ■ ■ — . . . . i . , , — I 
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FIG. 2: Relative differences between the 21cm power spectra for A-CDM and f(R) models, with fi = (i.e. for modes orthogonal 
to the line of sight). Power spectra are calculated at z = 11, assuming a neutral fraction xu — 0.9, accordingly to Ref. [60]. From 
top to bottom curves, the model parameter mo is respectively 5 x 10 _5 Mpc _1 (red), 10 _4 Mpc _1 (dark blue), 5 x 10 _4 Mpc _1 
(yellow), 10" 3 Mpc _1 (green) and 2 x 10" 3 Mpc _1 (blue). The grey band corresponds to the expected errors on the power 
spectrum measurements for the considered FFTT experiment. Errors are due to the cosmic variance at large scales and grow 
exponentially at small scales due to the angular resolution of the telescope. 



2. Symmetron 

The relative differences between the 21cm power spectra at z = 11 for the symmetron and the A-CDM models are 
plotted in Figs. 4 and 5, respectively for kn = and kj_ varying, and for k± = O.lMpc -1 and fc|| varying. At large 
scales, as for the f(R) model, one gets /3 2 /(l + m 2 a 2 /k 2 ) — > in Eqs. (47) and (48) and the symmetron cannot be 
distinguished from the A-CDM model. At small scales, /3 2 /(l + m 2 a 2 /k 2 ) — > f3 2 , inducing a scale invariant shift of the 
21cm power spectrum amplitude. 

In the case of orthogonal modes to the line of sight, Fig. 4 shows that the transition between these two regimes 
occurs in the range of observable scales for 10 _2 Mpc -1 < too < lOMpc -1 . The magnitude of the shift is controlled by 
/3q and by the redshift z„ from which modifications of gravity are triggered. For (3q ~ 1, modified gravity effects could 
be detected by the FFTT up to ~ 14, i.e. just before the reionisation. However, only values of /3q of the order of 
unity or higher will be detectable. 

When parallel modes are probed with k± ~ O.lMpc -1 and fcii < 2Mpc _1 , the transition regime is detectable up 
to mo ~ 200Mpc~ 1 (see Fig. 5). Moreover, the amplitute of the shift increases since the 21-power spectrum goes like 
(1 + /i 2 ) 2 with fi sa 1.1 at the smallest scales (it is larger than unity due to modified gravity effects in Eq. 11). As a 
result, the model signatures could be detectable up to z* ss 12, i.e. just before the redshit of observation, provided 
(3q w 1 or higher. 

Signatures on the matter power spectrum today can also be important for the symmetron model, as shown in Fig. 8 
for j3 = 0.5, too = O.lMpc -1 and z» = 20. Typically, values of m < lOMpc -1 with /3 — 1 and z* > 10 are already 
ruled out by observations. Nevertheless, as noticed above for the f(R) model, the matter power spectrum is limited 
to scales k < 0.3Mpc _1 and as a consequence its ability to probe large values of Too is reduced compared to the 21cm 
signal. 
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FIG. 3: Relative differences between the 21cm power spectra for A-CDM and f(R) models, with k±_ = O.lMpc -1 and k\\ varying 
up to 2Mpc _1 (cut-off introduced to avoid non linear effecs). As in Fig. 2, power spectra are calculated at z = 11, assuming 
a neutral fraction ih = 0.9. From top to bottom, mo values are respectively 2 x 10 _3 Mpc _1 (blue), 10 _2 Mpc _1 (red) and 
2 x 10 _2 Mpc _1 (yellow). The grey band corresponds to the expected errors on the power spectrum measurements for the 
considered FFTT experiment. 



Local constraints for the symmetron model are satisfied provided mo /Hq > 10 3 [58]. The matter power spectrum 
in the linear approximation gives a bound on mo for /3rj ~ 1 and z+ f=a 20 that is of the same order of magnitude. This 
bound could be improved by observations of the 21cm signal. 

To summarise, provided that symmetron effects arc triggered at rcdshifts larger than the rcdshifts of observation, 
the 21cm signal is found to be promising to put stringent constraints on the symmetron parameters, and especially the 
scalar field mass today, than the matter power spectrum and the local test of gravity (m < 200 Mpc^ 1 for z+ = 20 
and po ~ 0(1), i.e. approximatively three order of magnitudes better than local test constraints). Combining those 
signals and methods could be also a way to break the degeneracy between the model parameters (especially /3q and 
Z+) by probing different stages of the evolution of the matter perturbations and different environments. 

3. Dilaton 

For dilaton models, m(a) = vn§a~ 2 and /3(a) = /3oa 3 in the matter dominated era, with typically /3q of order unity. 
One thus gets ft 2 / (l+m 2 a 2 / k 2 ) = j3oa 3 /[I+ttiq/ (a 2 k 2 )] in Eqs. (47) and (48) for the evolution of matter perturbations. 
At redshifts z ~ 10, i.e. during the reionisation, the coupling to matter was therefore about a thousand times lower 
than today, so that the model signatures on the 21cm power spectrum are indistinguishable from the A-CDM case if 
we want the matter power spectrum today to remain under control. For large values of the scalar field mass (typically 
m > 0.1Mpc~ 3 ), the difference with the A-CDM model is even more suppressed since rriQ/{a 2 k 2 ) 1. 

4- Chameleon models 

The dilaton and f(R) models are particular cases of generalised chameleon models, for which /?o, mo, r and s can 
vary. In Figs. G and 7, we have calculated the 21cm power spectrum at reionisation for various values of /3q, r and s. 
Increasing /3q or s increases the relative difference with the A-CDM model while increasing r implies a reduction of 
this difference. 

In the case r > 3, s < and /3o ~ 0(1), the best constraints from local tests of gravity come from the galaxies, 
imposing ttiq/Hq > 10 3 , i.e. approximatively mo > O.lMpc -1 . This bound increases by several orders of magnitude 
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FIG. 4: Relative differences between the 21cm power spectra for A-CDM and symmetron models. As in Fig. 2, the power 
spectra are calculated at z = 11, assuming xu = 0.9, for wavelength modes orthogonal to the line of sight. Plain curves are for 
z* — 20, fio = 1, with mo varying. From top left to bottom right, mo takes the values 10 _2 Mpc (red), O.lMpc" 1 (green), 
IMpc -1 (yellow) and 10Mpc~ (blue). The two dashed curves are for /3o = 1, mo = O.lMpc - with z* varying. The top yellow 
one is for z* = 17, the bottom green one is for z* = 14. The dotted curve is for z* = 20, mo = O.lMpc -1 , /3o = 0.5. The 
parameters fio and z* are observed to control the magnitude of the shift between symmetron and A-CDM at small scales, 
whereas the parameter mo controls the scale of the transition regime. The grey band corresponds to the expected errors on the 
power spectrum measurements for the considered FFTT experiment. 
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FIG. 5: Relative differences between the 21cm power spectra for A-CDM and symmetron models, at z = 11, assuming xn = 0.9, 
for wavelength modes with k± — 0.1 Mpc - . Plain curves are for z* = 20, p = 1 and m = 10/100/1000 Mpc" 1 (respectively 
the blue, red and yellow plain curves). The two dashed curves are for /3o = 1, mo = 10 Mpc" 1 , z* = 14 (top blue) and z* = 12 
(bottom green). The grey band corresponds to the expected errors on the power spectrum measurements for the considered 
FFTT experiment. 
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FIG. 6: Relative differences between the 21cm power spectra for A-CDM and chameleon models, with fc|| = 0. As in Fig. 2, the 
power spectra are calculated at z = 11, assuming ih = 0.9. Fiducial parameter values are mo = 

and s = 0. Plain curves are for s = 0.5 (top yellow), s = (blue) and s = —0.5 (bottom red). Dashed curves are for r — 3.5 
(bottom green) and r — 2.5 (top blue) and dotted curves are for fa = 2/V6 (top yellow) and p = 1/(2^6) (bottom red). The 
grey band corresponds to the expected errors on the power spectrum measurements for the considered FFTT experiment. 



if r < 3 due to the stringent constraints from laboratory experiments with cavities. If s > 1, solar system tests of 
gravity can give the best constraints with the condition given by Eq. (87). 

If we impose s = and r — 3 together with m a = O.lMpc -1 (i.e. about the limiting value imposed by galaxy 
constraints) , we have found that the coupling to matter today needs to be (3 > 20 to lead to detectable effects in the 
21cm power spectrum with fc|| = 0. This limit is lowered to (3q ~ 2 for modes with k± = O.lMpc -1 and ku varying. 
For negative values of s, the limit on (3q increases. If we take positive values of s and impose /3q — 1 as well as r = 3 
and take 



/O 1f)4s-4r+18 

mo = H f-^— (88) 

(i.e. the minimal value for the solar system constraints to be respected), the effects on the 21cm power spectrum are 
never observable. An identical conclusion can be drawn for the matter power spectrum. Increasing the value of r while 
the other parameters remain fixed reduces the relative difference with the A-CDM model. Taking r < 3 increases this 
difference but it is compensated by the fact that larger values of mo are required to respect the bounds from cavity 
experiments. 

To summarise, it appears that the case r = 3 and s = is the configuration for which the 21cm signal could reach 
the sensitivity of the local tests, and signatures could be detected with a coupling to matter /3q > 2. Varying r and 
s leads to more stringent constraints from local tests so that such configurations should be very difficult to probe 
efficiently with 21cm experiments if the coupling to matter today is of the order of unity or lower. 



C. Probing /? 7 via the variation of a 

Due to quantum effects such as the presence of heavy fermions, the scalar field <p can be coupled to photons 

Sgaugc = -^2 J ^X^Bp^F^F^ , (89) 

where g is the bare coupling constant and Bp((j)) = 1 + /? 7 ft4</> + We consider the coefficients f3 and /3 7 as free 

parameters, even if they can be related depending on the model. The coupling to the electromagnetic field leads to a 
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FIG. 7: Relative differences between the 21cm power spectra at z — 11, assuming ih = 0.9, for A-CDM and chameleon models, 
with k± = O.lMpc -1 . Model parameters for the plain curves are mo = O.lMpc -1 , r — 3, s = and /3o = 5 (top blue), /3o = 3 
(red) and /3o = 1 (bottom yellow). The dotted curve is for mo = O.lMpc -1 , s — 0, j3o — 5 and r — 3.2 . The dashed curve is 
for s — 0.5, r = 3, /3o = 5 and mo = 0.8Mpc -1 (from Eq. 88). The grey band corresponds to the expected errors on the power 
spectrum measurements for the considered FFTT experiment. 



time dependence of the fine structure constant a [57, 58], 



— = -B F {4>) , 



(90) 



so that if we assume that Bf(<P) ~ 1, we get 



(91) 



The 21cm radiation is very sensitive to the variation of the fine structure constant, as noticed in Ref. [74] 3 . Indeed, 
the brightness temperature depends on the Einstein coefficient Aiq, itself depending on a. More precisely, Aio = 
27rcv^2i^p/(3c 4 7fte)' and j/21 <x a 2 Roo oc a 4 , where Roo is the Rydberg constant. Therefore the 21cm brightness 
temperature Tb oc Aiq/v^i oc a 5 , and the 21cm power spectrum Pat b <x aW - As a consequence, the relative error on 
a from a given experiment will be ten times better than the relative error on the 21cm power spectrum, 



Ac* _ 1 APaTb 
a 10 Pat b 



(92) 



From this argument, one expects that measuring the 21cm signal could be an excellent way to probe the variation of 
the fine-structure constant, and by extension a probe of the coupling /3 7 . However, in the context of modified gravity 
and a signal from the reionisation, the situation is not so ideal. The coupling to photons not only affects the 21cm 
brightness temperature through the time variation of the fine structure constant but also through the time-dependent 
fermion masses. Moreover, it is not clear how variations of a and fcrmion masses can affect the reionisation process 
itself, and thus the time evolution of ih, that is fairly unknow even in absence of modified gravity effects. 

From these considerations, it appears nearly impossible to constrain the value of /3 7 from single redshift observations 
only. Nevertheless, in the context of 21cm tomography over a broad range of redshifts (typically from the early stages 



3 Notice however that the results of Ref. [74] has been disputed and have lead to a dispute on how to calculate consistently the effects of 
the time variation of a on the 21cm signal [75, 76]. 
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FIG. 8: Matter power spectrum in the linear approximation for the A-CDM model (blue curve), for the f(R) model with 
mo = 2 x lCP^Mpc -1 (red curve), for the symmetron model with z* = 20, /3o = 0.5 and mo = O.lMpc -1 (yellow curve) and 



for the dilaton model with mo = 10 Mpc and 



1 (green curve). SDSS data points and error bars are also plotted. 



of the reionisation to the period following its completion, 12 > z > 2), the evolution of the mean ionized fraction could 
be reconstructed and combining high and low rcdshift measurements could provide a natural way of constraining the 
variation of a. Below, we estimate a bound on /3 7 that could be established from 21cm observations with the FFTT, 
under the assumption that the mass variation of fermions can be neglected in the brightness temperature. 

Assuming that the 21cm power spectrum amplitude will be measured with a percent level accuracy, one can constrain 
|Aa|/o; < 10~ 3 and a bound on the parameter /3 7 can be derived. From Eq. (91), one gets 



< 10" 



1 



(93) 



where A(f> is the scalar field variation during the redshift range of observation of the 21cm signal. It can be calculated 
from Eq. (70). For the f(R), dilaton and chameleon model, one gets 



9/3o#o 2 ^: 



2r- 



mg(2r 



2) 



2 ) 



(94) 



are respectively the maximal and minimal value of the scale factor in the range of the 21cm 
observations. For instance, for the redshift range given above and the f(R) model, one can obtain 



< 0.3^ 



(95) 



With nia/Ho ~ 10 , one gets /3 7 < 10 , which is an intermediate value between the present observational bound 
/3 7 o < 10 11 and the much tighter bound /3 7 o < 0.1 derived in [58] from the best experimental bound on the variation 
of a. The bound on /3 7 remains at a similar order of magnitude for dilaton and chameleon models, for reasonable 
values of r and s and as long as /3o = 0(1). 
On the other hand, for the symmetron model, 
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As an example, for 



20, one gets 



n 4 A(j) ~ 5 x 10 : 



(96) 



(97) 
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and thus 



|/3 7 | < 2 x 10~ 7 



777 ; 







(98) 
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With mo/Ho « 10 3 and /3q « 1, one gets /3 7 < 0.2, which is much tighter than the bound obtained for f(R), dilaton 
and chameleon models. 



The 21 cm line can in principle be used to probe the evolution of the matter perturbations over a wide range of 
redshifts, typically from the dark ages up to the completion of the reionisation. Observing the 21cm cosmological signal 
should therefore further our understanding of the evolution of the Universe and it is thus important to investigate 
the predictions of different cosmological models on the 21cm thrce-dimcntional power spectrum. In this paper we 
have considered modified gravity models with a screening mechanism and study their signatures on the 21cm power 
spectrum at rcionization. 

Our archetypical experiment is the Fast Fourier Transform radio- Telescope (FFTT) , consisting of a one kilometer 
side square of dipole antennas, that is designed especially for the detection of the 21cm power spectrum at reionisation. 
While the current and next generation of giant radio-telescope are expected to have a limited interest for cosmology, 
the ability of a FFTT-type experiment to put strong constraints on the various cosmological parameters has been 
demonstrated in [60]. 

We have investigated modified gravity models with a screening mechanism using a unified paramctrisation whereby 
the models are characterised by the scale dependence of the coupling to matter and the mass of the scalar field. This 
has previously been shown to encapsulate all the effects of modified gravity with a screening mechanism [57] and has 
been used to investigate the effects of such models on large scale structure [58] . 

For a fixed redshift, the consequences of modified gravity are important on small scales, within the Compton 
wavelength of the scalar field mediating the deviation from General Relativity. At the time of reionisation, such small 
scales could be in the linear regime of perturbations, where modifications of gravity enhance the growth of structure 
in a maximal way, whereas they are in the non-linear regime at the higher redshifts of the large scale structure 
probes, where screening effects take place and therefore reduce the magnitude of modified gravity effects. Hence 21cm 
cosmology offers a prime possibility of observing modification of gravity unhampered by screening effects. However, 
it is important to notice that with our specifications of the FFTT, such small scales can be only probed through 
wavelength modes parallel to the line of sight. 

The parametrisation of [57] is used here in the context of generalised chameleon, dilaton and symmetron models, as 
well as f(R) gravity, which is a special chameleon case. For all these models (except for dilatons), we find that 21cm 
observations at reionisation should constrain them more tightly than the large scale structures and CMB observations. 
The 21cm signal appears to be also a good discriminator of modified gravity models. 

Some models are already tightly constrained by the tests of gravity in the solar system, in the laboratory and in 
galactic environments. If these constraints are imposed, predictions for 21 cm cosmology are very similar to that of 
ACDM and signatures of modified gravity should be hardly observable. However, in the case of generalised symmetron 
models, strong constraints could be established with 21cm observations. 

We have also considered the effect of the scalar field coupling to photons in modified gravity models. This coupling 
arises naturally from the conformal anomaly [77] and gives rise to a time variation of the fine-structure constant. Since 
the 21 cm line is very sensitive to any variation in the fine structure constant, it can be used to probe the coupling 
to photons. This has enabled us to forecast tight bounds on this coupling on the basis of a 21cm experiment covering 
redshifts going from the beginning to the completion of the rcionization. 

More accurate predictions will require sophisticated MMCM and Fisher matrix techniques. This will allow us to 
consider the whole ensemble of wavelength modes accessible to the experiment and to perform a multi- redshift analysis 
in order to set precise bounds on the various model parameters, taking account the possible degeneracies with other 
cosmological and reionisation parameters. This is left for future work. 
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